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Abstract. A new approach is suggested for the study of geometric symmetries in general 
relativity, leading to an invariant characterization of the evolutionary behaviour for a class of 
Spatially Homogeneous (SH) vacuum and orthogonal 7— law perfect fluid models. Exploiting the 
1+3 orthonormal frame formalism, we express the kinematical quantities of a generic symmetry 
using expansion-normalized variables. In this way, a specific symmetry assumption lead to 
geometric constraints that are combined with the associated intcgrability conditions, coming 
from the existence of the symmetry and the induced expansion-normalized form of the Einstein's 
Field Equations (EFE), to give a close set of compatibility equations. By specializing to the case 
of a Kinematic Conformal Symmetry (KCS), which is regarded as the direct generalization of 
the concept of self-similarity, we give the complete set of consistency equations for the whole 
SH dynamical state space. An interesting aspect of the analysis of the consistency equations is 
that, at least for class A models which are Locally Rotationally Symmetric or lying within the 
invariant subset satisfying N a a = 0, a proper KCS always exists and reduces to a self-similarity 
of the first or second kind at the asymptotic regimes, providing a way for the "geometrization" 
of the intermediate epoch of SH models. 
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1. Introduction 

Due to their appealing geometric, kinematic and dynamical structure, Spatially Homogeneous 
(SH or Bianchi) models have received considerable attention in the last 3 decades. Apart from the 
obvious gain of a direct generalization of the standard Friedmann-Lemaitre cosmological model, one 
of the main reasons for this interest is the fact that the Einstein's Field Equations (EFE) are reduced 
to a coupled system of ordinary differential equations. Then by introducing the orthonormal frame 
formalism and expansion-normalized variables, in order to scale away the influence of the expansion 
from the overall evolution of the corresponding models, one exploits methods from the theory of 
dynamical systems to examine their behaviour at early, late and intermediate periods of their 
history. This approach has shown the significant role of the transitively self-similar SH models 
since they represent the past and future (equilibrium) states for the majority of evolving vacuum 
and 7— law perfect fluid models PEl- 

However, one can supplement this discussion with a broader consideration of the issues 
concerning the evolutionary era of SH models. The equilibrium state, whenever it exists, of evolving 
SH models is geometrically described by a model that admits a Homothetic Vector Field (HVF or 
self-similarity of the first kind) or a Kinematic Self-Similarity (KSS or self-similarity of the second 
kind) X which is defined according to: 

£ x m q = Su a Cx_h ab = 2ah ab (1.1) 

where h ab — g ab + u a u b is the projection operator normal to the timelikc unit vector field 
u a u a — — 1 and a, 5 =const. essentially represent the (constant) time amplification and space 
dilation |H EI EI d 

Therefore it will be enlightening to complete the geometric picture and find a relevant way 
to invariantly (although not uniquely) characterize the intermediate behaviour of the associated 
models, or to put it equivalently, one is tempted to ask: what was the nature of the generator of the 
self-similarity (either of the first or second kind) at the time of the intermediate evolution i.e. what 
kind of symmetry invariantly describes the intermediate behaviour of evolving vacuum or perfect 
fluid SH models? The natural and intuitive answer (although using heuristic arguments) is that 
one should expect that this type of symmetry must involve a generalization of the self-similarity. 
Consequently, the symmetry will represent the smooth transition mechanism between the evolving 
and the equilibrium states of SH models i.e. its "asymptotic behaviour" (into the past and the 
future) will be one of the well known symmetries. 

In order to achieve the above goal one could consider and study general symmetries taking 
into account the associated local diffeomorphisms and the intrinsic geometric structure of the 
spacetime manifold 8 . However, from a dynamical point of view, it appears natural to pursue a 
different direction by making use of the full set of EFE in order to augment and enrich the study 
of geometric symmetries. Clearly, this approach will provide us the necessary set of compatibility 
equations which can be studied in each subclass of SH (or even less symmetric) models. Motivated 
by the above discussion and the great success of studying SH models using elements from the 
theory of dynamical systems, in the present article we propose an alternative technique to study 
geometric symmetries that meets the aforementioned scope. This method fully incorporates the 
expansion-normalized orthonormal frame approach and give a transparent picture of how a specific 
symmetry "assumption" can be consistently endowed within a class of models or will produce 
further constraints, thus losing certain features of the general case. An illustrative and interesting 
example to which the preceding discussion applies is a recently presented new type of symmetry 
namely the concept of Kinematic Conformal Symmetry (KCS). This type of symmetry is regarded 
as a consistent, with geometry, generalization of the KSS or the case of conformal transformations 
preserving at the same time the causal structure of the spacetime manifold [HlEIj. Therefore, as a 
first step towards an invariant and geometric description of SH models, it seems adequate to study 
the implications from the existence of a KCS. 

The structure of the paper is as follows: in section 2, after a brief review of the basic elements 
of the 1+3 orthonormal frame scheme and the introduction of expansion-normalized variables, 
we specialize our study to SH geometries and give the tetrad/expansion-normalizcd form of the 
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kinematical quantities that define a generic symmetry. The determination of the corresponding 
expressions for the case of a KCS is treated in section 3. Since any type of symmetry assumption 
lead to geometric constraints which pass, through the EFE, into dynamics, we also determine 
the dynamical restrictions as follow from the associated integrability conditions. This set of 
equations, together with the expansion-normalized form of the EFE, constitute a close system of 
compatibility equations which can be used to check the consistent existence of a KCS in all the SH 
cosmologies. In section 4, special attention is given to some models of class A in which we study the 
consistency of the system of equations. We find that, at least models which are Locally Rotationally 
Symmetric (LRS) or lying within the invariant subset satisfying N a a = 0, always admit a KCS. 
A side result is that, at the asymptotic regimes, the KCS reduces to HVF or a KSS which shows 
that, in principle, the existence of a four dimensional group of KCS invariantly characterizes the 
intermediate behaviour of the associated evolving vacuum and perfect fluid models. Finally in 
section 5 we give some concluding remarks and comments and we discuss possible extensions and 
further advances which can be made in the direction of a satisfactory and conclusive answer to the 
"geometrisation" of the evolutionary behaviour of SH (or less symmetric) models. 

Throughout this paper, the following conventions have been used: spatial frame indices are 
denoted by lower Greek letters a,f3,... — 1,2,3, lower Latin letters denote spacetime indices 
a, 6, ... = 0, 1, 2, 3 and we use geometrized units such that 8nG = c = 1. 



2. The generic symmetry in SH models 

2.1. The 1+3 orthonormal frame formalism 

Spatially Homogeneous models are specified in geometric terms by requiring the existence of a Q% 
Lie algebra of Killing Vector Fields (KVFs) X Q with three-dimensional spacelike orbits S. This 
implies the existence of a uniquely defined unit timelike vector field u a {u a u a = —1) normal to the 
spatial foliations S: 

1 9 

U[a;b] = = U a . b U b -Cufjab = Ua:b = " a fc + -h ab (2.1) 

where o- a t,,8 = u a .j,g ab : h a b — g a b + u a Ub are the kinematical quantities associated with the u a 
according to the standard 1+3 decomposition of an arbitrary timelike congruence |llj . Because 
u a is irrotational and geodesic, there exists a time function t(x a ) such that u a = 81 i.e. each value 
of t essentially represents the hypersurfaces S. 

As far as the dynamical structure is concerned, from a cosmological point of view, it is sufficient 
to focus our study on SH models with an orthogonal perfect fluid assuming a 7— law equation of 
state. Therefore the EFE become: 

37 — 2 2 — 7 

Rab = ^ P U a u b H —pKb (2.2) 

It follows that the timelike vector field u a is identified with the average fluid flow velocity. 

On using the orthonormal frame formalism in SH cosmologies the starting point is to introduce 
a set of four linearly independent vector fields {e a } and their dual 1-forms {w a } which are invariant 
under the three-parameter group of isometries: 

pC a ,e o ] = 0, [X a ,u> a } = 0. (2.3) 

Then, by performing a time-dependent rescaling of {e a }, we can write (locally) the metric tensor 
in a manifestly Minkowskian form: 

ds 2 =r lab iu a (t,x a )uj b (t,x a ). (2.4) 

In this case the invariant vector fields {e a } and the connection forms T a bc uj c (where V ec e;, = 
r a hc e a ) satisfy the commutation and the first Cartan structure equations: 

[e a ,e fc ]= 7 c ab (i)e c , dio a = T\ c u b A ^ = ~y a bt J' A " c (2.5) 
where d and A are the usual exterior derivative and exterior product of 1-forms respectively. 
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via: 

r abc = i [v adl d cb + i lbd i d ac - vc d i d ba ] ^i a bc = - (r a bc - r a cb ) . (2.6) 

Furthermore, the requirement of the constancy of the metric under covariant differentiation implies 
that, r (afc)c = where T abc = rj d a^ d bc - 

The above definitions suggest that the tetrad form of the covariant derivative of every tensor 
field is written in the well-known way 

V. A'"//.. = e c (K a z.) + + ... - T d ic K a ^ - ... (2.7) 

where e c (K a £ ) is regarded as the directional derivative of the functions K a - \ along the vector 
field e c . 

Because in SH (vacuum or perfect fluid) models there always exists the preferred and well 
defined timelike vector field u a associated with a congruence of curves normal to the three- 
dimensional spacelike orbits S of homogeneity, it is natural to select it as the timelike frame 
vector i.e. eo = u. It follows from equation l|2.5l) that the kinematical quantities of the timelike 
congruence u are directly related with the commutation functions J c ab according to [111 112) : 

J° Oa = u a = 0, 7 % = -2e Q/37 ^=0 (2.8) 

l\ a = -\^ a -^ a +^ oll ^ (2.9) 

where VL 1 is the local angular velocity of the spatial frame with respect to a Fermi-propagated 
frame along eo. On the other hand, the spatial components of 7 c afc are decomposed as: 

7°^ = a p 5 a 1 - a 7 S a p + e Mp n ap (2.10) 

leading to Bianchi class A and B models according to whether the quantity ap vanishes or not. 
Combining equations (|2.6() and 12.8() - (|2.10() we easily find: 

Q 

T^Oo = + a °<Pi ToqO = ^> r Q/ 3 = £ Q/ 3 7 f2 7 (2-11) 

r Q /3 7 = 2a[ a (5 / 3] 7 + e 7p [ a n p /3 ] + -e a f} p n p 1 (2-12) 

The complete set of the gravitational field equations is expressed in terms of the shear components 
cr Q( a, the expansion 9 and the spatial curvature quantities ap, n a p, by utilizing the Ricci identity 
for u a , the Jacobi identities for the frame vector fields and the Bianchi identities. 

2.2. Expansion normalized variables 

Of particular importance in the exploration of the asymptotic dynamics of SH models, is the 
reformulation of the complete set of orthonormal frame equations, as an autonomous system of first 
order ordinary differential equations. This can be done by defining a set of expansion-normalized 
(dimensionlcss) variables which results the decoupling of the evolution equation of H = 8/3 from 
the rest of the evolution equations: 

dt 1 dH . . . „ . 

Tt = h> ^ = - (1 + < ^ (2 - 13) 

where q, H are the deceleration and Hubble parameter respectively and r is the dimensionlcss time 
variable. 

The complete set of equations can be written in the form |13| : 

J? a p = -(2-q) S Q/3 + 2e' lv { Yl^R v - S a p (2.14) 
Kp - QN afi + 2T,\ a N 0)lx + 2e p \ a Np )p R u (2.15) 
A' a = qA a -W a A fl + e^ a A fl R l/ (2.16) 
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Q' = Q [2q - (37 - 2)] 



(2.17) 



where: 




(2.20) 



(2.18) 



(2.19) 



where: 



K=\- (2N\N Pl - N\N a p) 6^ + A 1 A 1 



(2.21) 



and the deceleration parameter is given by the relation: 



g=-S^E^ + -fi[(3 7 -2)] 
= 2(l- J fT) + ifi[3 ( 7 -2)]. 



(2.22) 



2.3. The generic symmetry in expansion-normalized variables 

The folklore for investigating the implications of the existence of geometric symmetries in General 
Relativity can be divided into two main categories. The first category is a geometric approach in 
which we study symmetries taking into account the holonomy group structure of the spacetime 
manifold together with the associated local diffeomorphisms 8 . In the second category one 
formulates the necessary and sufficient conditions, coming from the existence of the symmetry, in a 
covariant way and study their consequences in the kinematics and dynamics of the corresponding 
model ^31 • Of course one could also deal directly with the resulting system of partial differential 
equations (pdes), presupposing a specific geometrical and dynamical configuration which render the 
symmetry equations to be more tractable. Obviously this approach undergo many disadvantages 
and pathologies. One of the serious stumbling blocks is the fact that as the generality of a model 
is increased (i.e. the underlying geometric structure of the model is less symmetric than the SH 
geometry) the symmetry pdes are progressively non-linear and very often lead to solutions of the 
EFE which are immediately ruled out physically. 

Here we suggest an alternative approach for the study of geometric symmetries which 
fully exploits the well-established orthonormal frame formalism in terms of expansion-normalized 
variables. Although this technique will be applied to a specific symmetry "assumption" (as we 
shall see in the next section this is not really an assumption, at least for a class of models, but a 
consequence of the complete set of EFE) it can be used in a more general context for the study of 
other types of important symmetries |15| . 

Let us consider an arbitrary vector field X and express its components in terms of the frame 
vector fields e a : 



where A, X a are continuously differential functions of the spacetime manifold. 

The first derivatives of X are decomposed into irreducible symmetry kinematical parts 
{ip, H ab , F ab } in the standard way: 



X = X a e a = Ae + X a e, 



(2.23) 



V b X a = ipg ab + H ab 



+ F ba 



(2.24) 



where 



4V; = W k X k , 



(2.25) 
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H, 



ah 



V fc X fc ) g ab 



(2.26) 
(2.27) 



are the conformal factor, the traceless symmetric part and the antisymmetric part respectively. 

Using the definition (j2~7j) and equations (|!HT|) . (|2~T2|) . (|2~23|) - (l2~T7|) . we find the tetrad 
analogue of the kinematical quantities. If we further invoke the expansion- normalized differential 
operators d a = e a /H in the general expressions we finally obtain: 



4V; = H [d (A) + d a (X a ) + 3A - 2A a X a ] 
H 



[-39 (A) + d Q (X Q ) + 3A - 2A a X a ] 



2H 0a = H [-d a (X) + 8 (X a ) - {8 ai + S Q7 ) X 1 + e^WXP] 



H a j3 



[Aj6 a (; 



H{d (p X a) + (5 a p + E Q/3 ) A 

2F 0a = H [<9 Q (A) + d (X a ) + (6 aj + S Q7 ) X 

1 



A( a <5 /3 ) 7 ] X 1 



e^WX 13 ] 



a/3 



H 



Ai„5, 



0h 



X' 



(2.28) 
(2.29) 
(2.30) 

(2.31) 
(2.32) 

(2.33) 



Equations (|2.28|I - H2.33|I represent the symmetry kinematical quantities in terms of expansion- 
normalized variables in SH geometries and can be used in order to have a first hint of how the 
dynamics affects on the geometry of SH models (or vice- versa). We note that, one could choose 
to define expansion-normalized symmetry kinematical quantities in a similar way we do for the 
shear and spatial curvature variables. This would be convenient for high symmetries where first 
or second derivatives of {ip, H a b, F a b} are involved. However for the cases we are interested in the 
present article the use of (|2.28|I - H2.33|I will be satisfactory. 

We conclude this section by pointing out that, because any type of symmetry assumption is 
described in terms of geometric constraints, these are inherited by the dynamics through the EFE 
()2.2|) . Therefore in order to visualize how the symmetry further interacts with the dynamics, it 
is necessary to determine the effect of the former on the Ricci tensor. By using the well-known 
commutation relation between the connection and the Lie derivative j!6j : 

= \ 9 ar [V c (£ X <?br) + V 6 (£x5cr) - V r (£x<?6c)] 



£xRab — V c 



(06) 



v (6 [r c Q) 



and the defining equation l|2.24|l . we can show after a straightforward calculation that: 

CxRab = -2V b V a ?A - 5afcV c V c ?A + 2V fc V (b if fc a) - V c V c if afc . (2.34) 

Essentially, the last equation represents a set of integrability conditions which can be used to check 
the consistent existence of every type of symmetry assumption. 



3. Generalized Conformal Symmetries in SH models 

Although there exists (and can be defined) a sufficiently large number of symmetries, the most 
important type of them (up to date) appears to concern the constant scale invariance of the metric 
represented by the existence of a proper HVF E| • F° r SH vacuum and perfect fluid models 
this is indeed the case due to the profound relevance of homothetic models with the equilibrium 
points of the SH state space |18M19ll2*Ull21| . 

Recently, a new type of symmetry has been suggested, the so-called bi-conformal 
transformations which can be interpreted as generalizing the concepts of the self-similarity and 
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the conformal motions. In the present work we will concern with an interesting subcase, that is, 
the so-called Kinematic Conformal Symmetry (KCS). In particular a smooth vector field X is the 
generator of a KCS iff the following relations hold 9 : 

£xu a = 8u a , Cxhab = 2ah ab (3.1) 

or, in terms of the metric: 

£ X ffafc = 2a# ab + 2 (a - S) u a u b (3.2) 

where a, 6 are smooth functions that we shall call symmetry scales and h ab = g ab + u a Ub is the 
projection operator perpendicular to the timelike congruence u a . 

Combining equations (|2.24l) and (|3.2I) we express the symmetry kinematical parts in the form: 

, 3a + 5 a — S , . , oN 

V = — ^ — , H ab = - (g ab + Au a u b ) . (3.3) 

It can be easily observed that when a — 5 the KCS reduces to a Conformal Vector Field (CVF) 
which, due to the equation (|3.1II . is necessarily inheriting i.e. the integral curves of u a are mapped 
conformally by the CVF X As a result the Lie algebra X of inheriting CVFs is always a 

subalgebra of the Lie algebra of KCS which shall be denoted as B. Moreover when the symmetry 
scales a, 8 are both (different) constants we recover the case of a Kinematic Self-Similarity 



3.1. Symmetry constraints 

Clearly there exists a direct dependence between the existence of a KCS, as well as any other type 
of symmetry assumption, and a specific cosmological model. This mutual influence is reflected in 
the induced geometric, kinematic and dynamic constraints which are imposed due to the intrinsic 
nature of the symmetry and/or as a consequence of the specific geometric and dynamical structure 
of the spacetime manifold. In the case of a KCS these constraints are derived from the general 



relations H2.28fl - I|2.31[l an( i the symmetry assumptions l|3.3|) : 

3a = H [d a {X a ) + 3A - 2A a X a ] (3.4) 

[-AHdo(X) + 3a + 5] =3 (a - 6) 5 = Hd (X) (3.5) 

- d a (A) + d (X a ) - (<5 Q7 + S Q7 ) A 7 + e aM ICX fi = (3.6) 

H{d((3 x a) + (<W + ^W) A - [A 7 5 al 3 - A( a 5p-)y] X 1 

+ e l5[a N 5 p) X^} ~ a5 aP = 0. (3.7) 



We should emphasize that the above set of constraints must be augmented with the associated 
Jacobi identities satisfied by the generators {X Q , X} that constitute the Lie Algebra B of KCS. This 
will require the determination of the dimension of B and the assumption that the former is finite 
since there are cases for which B is infinite dimensional [Sj- Nevertheless we shall not pursue the 
problem in full generality and we will restrict our considerations to the case where the dimension 
of B is finite and equal to four (together with the G3 of KVFs which can be seen as "trivial" 
KCS). In fact this assumption seems reasonable and not even restrictive, from a geometrical point 
of view, because we intend to describe the "intermediate behaviour" of a (proper) self-similarity 
of the first or second kind in a way that is identified with the presence of a KCS. In this case it 
can be shown from Jacobi identities 233 that the Lie bracket of a KCS with the KVFs is always a 
linear combination of the later. It turns out that the scalar A and the symmetry scales a, S (due 
to equation (|3.1J) ) are spatially homogeneous: 

d a X = d a a = d a S = 0. (3.8) 

In addition we point out that equation (|3.2|l is a consequence of (|3.1|l since the former suppresses the 
scale amplification of both time and space which is essentially represented by the latter. However 
the spatial homogeneity condition (|3.8|) ensures that the definitions (|3.2|) and l|3.1|l are equivalent. 
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3.2. Integrability conditions 

Due to the geometric character of the KCS to preserve the causal structure of the spacetime 
manifold, it is natural to expect that the associated transformation group will respect, to some 
level, the intrinsic properties of SH models. Therefore one should expect that the presence of a 
KCS will induce weaker constraints than the case of CVFs or the KSS. In order to confirm this, it 
will be required to determine the integrability conditions adapted to the case of SH models filled, 
in general, with a 7— law non-tilted perfect fluid. 

As we have shown, the scale functions and the trace i/j are both spatially homogeneous which 
implies that equation (|2.34l) . due to (|3.3() . can be written: 



£xRab — 3 



H 5 - 2a 



u a u b 



2{a-S) (3H 2 + H 



H 6a - 5 



h a b 



{ia - S\ a ab + {a - 5) (6Ha ab + 2a ab ) 



(3.9) 



where a dot "•" denotes differentiation w.r.t. u a . 

On the other hand the EFE (|2.2jl and the relations (|3.1|) give: 



£x.Rab 



3 7 



[X (p) + 2Sp] u a u b 



7 



[X (p) + 2ap\ h ab . 



(3.10) 



2 " 2 

The complete set of integrability conditions, coming from the existence of a KCS, follow by equating 

(E3 and Ennt : 

37-2 



H 16 - 2d 



2 

a + H 



X(p) + 26p] 



2(a~6) (SH 2 + +& + H Ua-S) 
(3a - S^j a ab + {a - 5) (6Ha ab + 2a ab ) - 



7 



X (p) + 2ap] 



0. 



(3.11) 
(3.12) 
(3.13) 



We observe that in the case of a CVF where a = 5, the last equation implies the well-known result 
ctGab = i.e. either the (necessary inheriting) CVF reduces to a HVF or the spacetime is the 
Robertson- Walker spacetime |24j . 

The system of equations (I3.11|I - (I3.13|) can be conveniently reformulated in expansion- 
normalized variables in order to append them in the autonomous set ()2.14|l - (|2.17|l and the symmetry 
equations (|3.4|I - H3.7|I . For simplification purposes we define the dimensionless symmetry scale 
functions: 



(3.14) 



Then, taking into account equations l|2.14|l and l|2.17|l . an appropriate combination of l|3.11|l - H3.13|l 
eliminates the second order time-derivatives and give evolution equations for a and S: 



4 [5' - (g + 1) 5] 



2 (a - «) (2 - q) 
(37 - 2) 5 + 3 (2 - 7) a - 6A7 



n 



(3a- 5\ - (g + 1) (sa-S) £ a/3 - 2 (a - 6^ S, 



a/3 



(3.15) 
(3.16) 



It is interesting to note that equation H3.16(l excludes the existence of a proper KSS in Bianchi 
vacuum or perfect fluid cosmologies with S a /3 7^ 0. The constraint S a p — is identically satisfied 
in Kasner type I models which is well known to admit a (proper) self-similarity of the second kind 
[SJ • This result implies that a KSS fails to be considered as a generic candidate to describe 
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the intermediate behaviour of general SH models, but rather one should explore the possibility 
for a symmetry, representing a direct generalization of the conformal motions, in such a way that 
its asymptotic behaviour is the self-similarity of the first or second kind. We shall demonstrate 
in the next section that the case of a KCS provides an evidence towards a satisfactory (but not 
conclusive) answer to this question, for a significant subclass of evolving SH vacuum and perfect 
fluid models. 



4. Application to SH models of class A 

Spatially homogeneous cosmologies of class A are defined by the condition A a = 0, which due 
to equation H2.20JI . implies that the shear and spatial curvature matrices E a p>N a p commute. 
Therefore T, a p,N a p have a common eigenframe and we can write: 

£ Q/3 = diag(£ 11 ,S 22 ,£ 22 ), N af3 = ding (Nu, N 2 2,N 2 2). (4.1) 

Furthermore the evolution equations H2.14J1 . I|2.15(l and equation j|4.1|l show that R a = which 
means that the common eigenframe of S a /3, N a /3 is also Fermi-propagated. Under these conditions 
the study of the class A models is considerably simplified and permits one to investigate the set of 
consistency equations constituting of the EFE l|2.14|) - (|2.17l) . the geometric constraints (I3.4I) - (|3.7|I 
and the integrability conditions (|3.15|l - (|3.16|l in a straightforward way. To illustrate the method 
that could be used for the consistency checking we outline some applications by restricting our 
study to type I, II and VIo models. 

However, before we proceed, it will be useful to give the corresponding analysis for the flat 
Friedmann-Lemaitre model T since it represents the past or future attractor for several SH models 
of class A. 



4--1. Flat Friedmann-Lemaitre j—law perfect fluid models 

It is well known that the (flat) Robertson- Walker spacctimc: 

ds 2 = -dt 2 + S 2 {t) (dx 2 + dy 2 + dz 2 ) (4.2) 

has a variety of ways for an invariant characterization of its structure. For example, kinematically, is 
defined by the vanishing of the shear, vorticity and acceleration of the preferred timelike congruence 
u a which, due to equations (|2.19(l and 12.22JI . implies: 

Sl = l, ?=^rr^ (4-3) 

On the other hand one can use geometric terms and describe the Friedmann-Lemaitre model by 
the existence of nine proper CVFs, one of which is parallel to u a |5Sl|27j. Although in the case of 
a 7— law perfect fluid model a proper HVF and KSS exists |3] (lowering the dimension of the Lie 
algebra of conformal motions to eight) this does not mean that we have exhausted all the possible 
(geometric) ways for the description of Friedmann-Lemaitre models. As a result one should expect 
that a KCS will also exists without imposing extra geometrical or dynamical restrictions (which is 
often the case for other types of symmetries). Indeed from the constraints ()3.4|1 - H3.7|) we find: 

3(5 -A) = d a (X a ), 5 = A' (4.4) 

d o {X a )-X a = 0, d (p X a) + (A - a) 5 a p = 0. (4.5) 

Using the set of equations l)4.3[l - l)4.5[l and the integrability condition (|3.15(l we can show that a 
KCS always exists in a Friedmann-Lemaitre model and the symmetry scale a is given in terms of 
the temporal component: 

(X-a)' = (q + l)(X-a) = ^(X-a). (4.6) 

From the above equations we determine the exact form of the KCS in the Robertson- Walker 
spacetimc: 

X = X(t)d t + c (xd x + yd y + zd z ) (4.7) 
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where: 

ait) = m+^lt (4 8) 

and the scale factor is given by S(t) = t 2 l Zl . We note that, after a suitable change of the basis of 
the KCS Lie algebra, we can set the constant c = which implies that X is also parallel to the 
fluid velocity u a . 

4-2. Type I models 

In this case N a p — — Sap and equation (|3.16|) gives: 

3a-<5 = S= ^^3a-<5 = c (4.9) 

where c is an arbitrary constant. 

In addition, equation l|3.5[) implies that: 

6 = A'. (4.10) 

We should pointed out that the symmetry constraints are necessary to ensure the existence of a 
KCS. However this does not imply that they will be preserved along the integral curves of the 
timelike vector field eo. Therefore we must propagate equations l|3.4|l - H3.7(l in order to retain the 
existence of a KCS in every spacelike hypersurface S. After a short calculation and the use of the 
commutator relations (|2.9(l . we obtain: 

A' - [q+ 1) A = a - [q+ I) a (4.11) 



where 



[(<? — 2) A + A' — (g + 1) A] S Q/3 + A Q(3 = (4.12) 



A Q ^ = S 7 ( Q A 7 ^ - S 7 ( Q X /3 )' 7 (4-13) 



and A 7 p = dfjX^ . 

From the act— component of l|4.12[l it follows that: 

( (? -2)A + A'-( (? +l)A = 0^A' = 3A (4.14) 

while the equation K a p = expresses the spatial components of the KCS in terms of the shear 
variables. 

In summary we have shown that every type I 7— law perfect fluid model always admits 
a four dimensional group of Kinematic Conformal Symmetries. An interesting feature of this 
result concerns the past "asymptotic behaviour" of the KCS. In particular, from equation 1|4.14[) 
we observe that A' — (<j"+l)A = (2 — q) A hence at the equilibrium point (q = 2) we have 
A' — (q + 1) A = i.e. a, S =const. and the KCS reduces to a proper KSS with A = const, x H . 
The future state of the KCS is treated similarly. The type I models, approach at late times the 
Friedmann-Lemaitre model so equation (|4.ti|l is trivially satisfied and the temporal component 
is given in equation l|4.14|l . This implies that a' — (q + 1) a 7^ a, 5 ^const. as expected. 
Therefore, one could argue that the existence of a proper KCS describes, in a geometric fashion, 
the intermediate behaviour of the evolving type I models. 

We conclude the type I case by giving the exact form of the KCS in local coordinates. The 
general 7— law perfect fluid solution can be conveniently written in the form yQ: 

ds 2 = _ A 2h~D dt 2 + J2 t 2p« A W/3-P<*) ( dx «f (4.15) 
a 

where the function A 2 1 = k + m 2 t 2 1 and k, m are constants. 
It follows that the (local) coordinate form of the KCS is: 

X = \{t)d t + (ck + c 3 ) xd x + c 2 yd y + c 3 zd z (4.16) 
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X(t) 



tkcA 



(p 2 + 2p 3 -l)(A-iA ) t) 
and the constant c 2 is given by: 

C3 (j>2 + 2>P3 - 1) + fee (p 3 - p 2 ) 

C2 = 

P2 + 2p 3 - 1 

where we have used the well-known relations: 
Pi +P2+P3 = 1, 



pI+pI+pI = 1 



(4.17) 



(4.18) 



(4.19) 



satisfied by the Kasner exponents p a . 

We also note that, at early times X(t) sw t and a, 5 =const. whereas at late times, k = and 
a, 8 ^const. which confirm the reduction of the KCS to a proper KSS and to the KCS (|4.7|l 
respectively. 

4-3. Type II models 

The Bianchi type II invariant set is characterized by the conditions Ni > and N2 — N3 = 0. We 
find convenient to collect the consistency equations as follow from (|2.14() - (|2.17() . (|3.4() - (|3.7|l and 



X' -(g + l)A = a'-(g + l)a 

= (A' - 3A) E af} + A afj + [2EhA-t + JVie 7 i (Q $ 



- \8l5l - \5l5l ) A -f\ 



4[5 / -(g + l)5] 



3a- £ 

1 



2 (a - 5) (2 - g) 
(3 7 - 2) 5 + 3 (2 - 7) a - 6A7 



(i + q ) (3a -s 



2 X 2 



X' 



12 



a u 



n = o, 



2XH 



3 7 -2 



O 



= (g - 2)£ a/3 - - \S 2 A - JV? 



N[ = (g + 2E 11 )JV 1 . 



(4.20) 
(4.21) 

(4.22) 

(4.23) 
(4.24) 
(4.25) 

(4.26) 
(4.27) 



Note that, in complete analogy with the type I models, equations (|4.2()|) and (|4.21(l are the result 
of the propagation of the symmetry constraints (|3.4|) and l|3.7|) along do. 

A quick observation can be made, due to the form of (|4.21|) or (|4.23|l . For example, the first 
equation implies^: 

(A' - 3A) S 22 = (A' - 3A) £33 = -~iV x 2 A (4.28) 
where for E 22 = the KCS turns to the KVF of Bianchi type II models. 

I We recall the traceless property of T, a a i.e. Sn + S22 + S33 = 0. 
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Equation Ij4.28(l means that the existence of a KCS is compatible only with type II models 
which are LRS i.e. only within the invariant subset S^(II). In order to determine the symmetry 
scale a we make use of the a/?— components of equation (|4.21|) . Then, from the spatial commutators 
(|2.1UI) we get X% t i — -^3,1 = and the consistency of the remaining set of equations is assured 
provided that: 

a = X (4£ 22 + 1) . (4.29) 

We have proved that every type II evolving vacuum or j—law perfect fluid model that belongs to 
the invariant subset (II) can be invariantly characterized by the existence of a four dimensional 
group of Kinematic Conformal Symmetries. 

Regarding the "asymptotic behaviour" of the KCS in type II models, it is straightforward to 
show that in S^(II)\q-q the following relations hold: 

6t — (q + 1) a oc — 1 
3a + 6\ - (q+ 1) (sa + Sj oc £| 2 - 1. 

Therefore, at the equilibrium points £22 = ±1 and iVi = 0, the KCS becomes a proper KSS as 
expected, since at the asymptotic regimes, all models within (II)\n=o approach some vacuum 
Kasner model |28| . 

On the other hand, it is well known that non- vacuum models S^~(7/)|n>o are all future 
asymptotic to the homothetic Collins model (II) for which we have proved that does not admit 
a proper KSS, and past asymptotic to a Kasner model /C or to the Fricdmann-Lcmaitrc model T . 
Consequently, every orbit in Si(II)\q>q joins two (first and second kind) self-similar equilibrium 
points, hence we expect that the KCS will be reduced to a proper HVF and a KSS, except to 
the case where the orbit lies in the one-dimensional unstable manifold of T . This model has past 
attractor the point T and the KCS will reduce to the associated KCS of the Robertson- Walker 
spacetime (equation l|4.7p ). 

Indeed, using equations (|4.24() - (|4.29|l we find: 

- b = A [Nf + 6S 22 (2S 22 - 1)] 
3S 22 

From the last equation we can show easily that — <5^j |p+(jj) = which implies that 
lim (a — S) =0. Similarly we can show that, at the equilibrium point JC the symmetry scales 

a, (5 become constants with a\fc 7^ 8\k. and at T we have a = XH, X' = 3A. 
The exact form of the KCS in local coordinates is found to be: 

X = X(t)d t + (lex + by) d x + cyd y + (cz - 1) d z (4.30) 

where c is constant and: 

The LRS type II metric is: 

ds 2 = -A 2 dt 2 + B(dx + bzdyf + C (dy 2 + dz 2 ) (4.32) 

for smooth functions B(t),C(t) of their argument. 

For example, in the LRS vacuum model (a special case of the general solution found by Taub 
the metric is Q]: 

ds 2 = -A 2 dt 2 + t 2pi A- 2 (dx + 4 Pl bzdy) 2 + t 2p3 A 2 (dy 2 + dz 2 ) (4.33) 

where A = (l + bH^ 1 ) 1 ' 2 and p a satisfy PJSJ) . 

The temporal component of the KCS is given by: 

t A 

X(t) = — 7 — - (4.34) 

w 2tA it - ( Pl - P3 ) A v ; 
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and reduces to a KSS for small and large values of the time coordinate (the KSS in plane symmetric 
Bianchi models has also been found in |3"H]). 

As a final remark we note that by defining the new time coordinate t — to = C(t)/ B(t), the 
KCS becomes X = 2ctdf + (2cx + by) d x + cyd y + (cz — 1) d z . Then at t — we have a a b = and 
equation 13.2fl implies c = i.e. the KCS is reduced to the KVF of type II models. 



4-4- Type VIq models 

Let us consider now the Bianchi type VIo invariant set which is defined by Ni = and -ZV2-ZV3 < 0. 
Propagating equation l|3.7[) and using (|2.9() we obtain: 



= (A' - 3A) S Q/3 + K afi + [2£ 22 jr + E\A e ] N 2 e l2(ct 6% 

+ [2S 33 X^ + E\A e ] JV 3 e 73 (a$) + (4.35) 

where S a p and are given in 12.18fl and l|4.13(l respectively. 
The a/3— components of l|4.35(l . after a short calculation, give: 

(E 22 - E 33 ) [28 [3 X 2] + X 1 {N 2 + N 3 )] = (4.36) 

(2E 22 + S 33 ) (2d [2 X n + X 3 N 2 ) = = (2E 33 + E 22 ) {-2d [3 X 1} + X 2 N 3 ) (4.37) 

which implies that S 22 = S 33 since it can be shown that for E 22 ^ E 33 the KCS is becoming a 
proper HVF. 

In this case N 3 = -N 2 and the system of equations (|2~TI|l - (|2~T7|l . EFty - $£7\ and ^J^-^HM 

is identically satisfied, provided that: 

<5 = A(1-2E 22 ) (4.38) 
and the function A is given in terms of the shear and spatial curvature variables (S 22 =/= 0): 

V = M^Hl. (4.39) 

OZj 2 2 

Note that for E 22 = the above system of equations implies A = = a = 8 and the KCS reduces 
to an isometry. Moreover, it is not difficult to show that: 

(a-8)\^ = -*^±» (4.40) 
(~ x\t 2A [A| + 3E 22 (E 22 + 1)] 

(a-<5j|n>o = [ -• (4.41) 

Due to equations (|4.4(J|) and (|4.41|) . the analysis of the asymptotic behaviour of the KCS is 
straightforward. The past and future attractors of the vacuum invariant subset are the equilibrium 
points E 22 = 1,N 2 = (the arc K-i) and E 22 = —1,N 2 = (the Taub point %) respectively 
and the KCS reduces to a proper KSS (with 3a + 5 = 4ip = 0) and a proper HVF. Regarding 
the non-vacuum models it has been shown that the future attractor is the Collins homothetic 
model Pi {VIq) for which (a — S)\ p +^ VIo ^ — and the KCS reduces to a proper HVF. The past 

equilibrium state of the invariant set S^(VIo) is either the LRS Kasner point Qi where the KCS 
becomes a proper KSS or the Friedmann-Lemaitre model T and the KCS is given in equation (|4.7|l 
with A' = 3A. 

These results suggest that, every type VIq evolving vacuum or 7— law perfect fluid model lying 
in the invariant subset S{(VIq) always admits a four dimensional group of KCS that reduces, at 
the asymptotic regimes, to a self- similarity group of the first or second kind except from a set of 
measure zero for which the KCS preserves its nature. 

For completeness we also give the form of the KCS for the general vacuum solution satisfying 
N a a = (an arbitrary multiplication constant has been omitted) |31l I32| : 

ds 2 = r x ^e e (-dt 2 + dx 2 ) + t {e 2x dy 2 + e- 2x dz 2 ) (4.42) 
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d t + c x d x + c 2 yd y + (2ci + c 2 ) zd z 



(4.43) 



where the symmetry scales are: 



(ci + c 2 )(4t 2 -l) 
(4t 2 - 3) 



5 = 



(ci + c 2 ) (4t 2 - 9) (At 2 + 1) 
(4i 2 - 3) 2 



(4.44) 



We can verify that, at early and late times, we have 3a+i5 ss and a— S ~ 0, so the expected limiting 
cases of the KCS is a self-similarity of the second or first kind§. Also to be noted is the apparent 
"singular" behaviour (due to the specific time gauge) of the KCS as t — » V3/2 (=> a a b — > 0). 
However, similarly to the previous case, we can show that at this value, the temporal component 
of X vanishes and the symmetry equations imply that c\ + c 2 = i.e. the KCS X — > d x — yd y + zd z 
reduces to the KVF of Bianchi type VIo models. 

5. Discussion 

As we have mentioned, it is customary to study a specific symmetry assumption from a geometrical 
point of view without taking into consideration the kinematical and dynamical structure of the 
corresponding model. As a consequence, the effects on the dynamics of the symmetry constraints 
are hidden and usually produce unphysical results. Up to date, self-similarity of the first kind 
appears to be the only symmetry condition with a transparent physical nature since it represents 
a geometrization of the asymptotic (equilibrium) state of general models. Therefore, in order to 
complete the dynamical picture, it was of interest to seek and find a symmetry that could be used 
effectively as a consistent tool for the invariant description of the intermediate behaviour of general 
models. 

In the present article we have proposed a new technique of studying geometric symmetries by 
fully exploiting the 1+3 orthonormal frame scheme and the introduction of expansion-normalized 
variables. The analysis of the complete set of consistency equations IJ2.14Jl - IJ2.17Jl . J3.4Jl ~ H3.7JI and 
(J3.15Jl - (J3.1fiJl has revealed a novel feature of a large class of SH models that is summarized in the 
following: 

Proposition 1 Evolving SH models of Bianchi type I, LRS models of Bianchi type II and type VIq 
models within the associated invariant subset N a a = 0, are geometrically described by the existence 
of a four dimensional group of KCS that is reduced to a self- similarity group of the first or second 
kind at the asymptotic regimes, except from a set of measure zero for which the properness of the 
KCS is preserved. 

However, as we have seen, even in the exceptional cases the corresponding equilibrium state 
is the Friedmann-Lemaitre model T which we have proved that always admits a proper KCS, 
supplementing the geometric properties of the standard cosmological model. 

At first sight, the existence of a proper KCS in SH models appears somewhat surprising, at 
least from a dynamical point of view. This is mainly because the interaction mechanism between the 
geometric "assumption" of a KCS and the dynamical behaviour of SH models is not, conceptually, 
apparent. However a closer look on the structural properties of SH models indicates a possible 
qualitative interpretation of this interaction. In particular, the full set of non-linear EFE can 
be seen as the perturbed version of the associated linearization of equations IJ2.14JI - IJ2.17JI . at the 
vicinity of a hyperbolic equilibrium point. Accordingly we may interpret the generator of a KCS as 
representing the perturbation (to some order) of the corresponding generator of the self-similarity 
transformation group of the first or second kind. Eventually, this observation will enable us to 
geometrize the majority of the concepts and techniques that are used in the theory of dynamical 
systems with a view to optimize and efficiently elaborate the results from the qualitative study of 
general cosmological models. 

§ We recall here that the time coordinate t does not represent the proper clock time as t — > + or t — > +oo. 
However, the freedom of choosing the time gauge, ensures a similar asymptotic behaviour for the symmetry scales. 
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Although the existence of a KCS signifies the physical ground for a first promising attempt 
towards a "geometrisation" of the evolutionary behaviour of SH models, we do not allegate that 
the KCS (uniquely) characterizes the intermediate epoch of the totality of SH models. In fact, 
because a KCS possess two arbitrary spatially homogeneous functions (the symmetry scales), we 
expect that only SH models with two essential degrees of freedom will exhibit a proper KCS. This 
conjecture is confirmed by the Proposition 1 in which all the exact solutions found to admit a 
proper KCS belong to this class. The main reason, that provides an interpretation for the possible 
connection between exact solutions and the existence of KCS, appears to involve the so-called 
"hidden" symmetries of the SH cosmologies. Indeed, using the Hamilton- Jacobi reformulation of 
the EFE in which all the dynamical picture is encoded in one geometric object, it has been shown 
that all known exact solutions with two degrees of freedom are associated with a specific kind 
of "hidden" symmetry namely the existence of a Killing tensor symmetry of the Jacobi metric 
that generalizes the corresponding cyclic variables and the Hamilton- Jacobi separability [331 I34| . 
Therefore it will be of interest to extent the analysis to the rest of the Bianchi models in order to 
see if the existence of a KCS is a general feature of the two-dimensional SH invariant subset i.e. if 
it is related with the above type of "hidden" symmetry of the Bianchi cosmological models |15j . 

We should remark that, although we have mainly focused our study to Bianchi types with 
clear and simple past and future equilibrium states, one could apply the approach presented in 
this paper to models with more complicated dynamical structure e.g. models with oscillating or 
diverging asymptotic behaviour near to the past or future attractor. The non-vacuum Bianchi 
VIIo invariant subset provides an interesting example since it is well known that the associated 
kinematical variables are unbounded and do not approach any equilibrium point into the future i.e. 
non- vacuum Bianchi type VIIo models are not asymptotically self-similar As a consequence 

one should expect that a KCS does not exist in those models due to the non-existence of a self- 
similar model as future attractor. However, a preliminary analysis has shown that a KCS does 
exist in LRS type VIIo models which is never reduced to a HVF or a KSS suggesting that 
the concept of the KCS represents not simply a perturbed version of the self-similarity group but 
a generic property (in the spirit of [331 I34| ) of the two-dimensional invariant subset of the SH 
cosmological models. 

Clearly, the case of higher dimensional invariant subsets requires further investigation. 
Assuming the existence of several proper KCS, will not solve the problem for models with three or 
more essential degrees of freedom, since the condition (|3.1[l implies that the symmetry scales are 
always spatially homogeneous restricting the dimension of the Lie algebra of KCS in SH models 
to four. Therefore the question of determining the symmetry which invariantly describes the 
whole set of SH models is still open. Nevertheless, the implications of the above results, enforce 
the important role which may play a specific (still unknown) general symmetry as an effective 
geometric implement for the invariant description of general cosmological models and not only 
as a simplification rule towards the determination of exact solutions with ambiguous (or even 
without any) physical meaning. A closely related issue is how the constraints, coming from the 
presence of the general symmetry, could reveal a path of constructing the general (whenever is 
possible) solution of the corresponding cosmological model. We expect that, the approach of 
studying generic geometric symmetries in spacetime developed in this article, can be applied to 
more general geometric setups leading to a more efficient qualitative and analytical study of general 
vacuum and perfect fluid models. 
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